SIMULTANEOUS GENERATION FOR ZETA VALUES BY THE 

MARKOV- WZ METHOD 
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Abstract. By application of the Markov- WZ method, we prove a more general form 
of a bivariate generating function identity containing, as particular cases, Koecher's 
and Almkvist-Granville's Apery-like formulae for odd zeta values. As a consequence, 
we get a new identity producing Apery-like series for all C(2?i + 4m -I- 3), n, m > 0, 
convergent at the geometric rate with ratio 2"^". 



1. Introduction 
The Riemann zeta function is defined by the series 

oo _ 

^(^) = E;^' Re(3)>l. 

n=l 

Apery's irrationahty proof of ({3) [13] operates with the faster convergent series 

first obtained by A. A. Markov in 1890 [lOj. The general formula giving analogous series 
for all ({2s + 3), s > 0, was proved by Koecher [7] (and independently in an expanded 
form by Leshchiner [S]) 

(1) ^C(2s + 3)x2^ = ^^^^^--^ = -^-— 2^ n [^-^J- 

s — k — 1 k — 1 \ K / iTi — 1_ 

The similar identity generating fast convergent series for all C(4s + 3), s > 0, which for 
s > 1 are different from Koecher's result ([1]) was experimentally discovered in [3] and 
proved by G. Almkvist and A. Granville in [Ij 

s=0 k=l k=l \k) m=l ^ ^ 
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There exists a bivariate unifying formula for identities ([T]) and ([21) 

(3) 1^ f.A _ ^2f.2 _ yA - 2 1^ A;4_a,2p_^4 11 I ^4_a;2^2_^4 

fc=l ^ k=l \k) ^ m=l ^ ^ 

which was first conjectured by H. Cohen and then proved by D. Bradley [5J and, inde- 
pendently, by T. Rivoal [14j. This identity implies if y = 0, and gives ([2]) if x = 0. 
The proof of ([3]) relies on Borwein & Bradley's method [3] and consists of reduction 
of ([3]) to a finite non-trivial combinatorial identity which can be proved on the base of 
Almkvist & Granville's work [1]. 

In this short note, we prove a more general form of ([3]) by application of the Markov- 
WZ method (see [IHl ESI [HI E])- Let us notice that Koecher's identity ([T|) and similar 
Leschiner's and Bailey-Borwein-Bradley's identities [21 [1] generating rapidly convergent 
series for even zeta values ({2s + 2) can be proved with the help of the WZ method (see 
[6] for more details). 

Theorem 1. Let a,b be complex numbers, with \a\ < 1, |6| < 1. Then for arbitrary 
complex numbers Aq,Bo, Cq we have 



- (F - a2)(P - 62) ^ n:;=iK - - h^) 



with 



_ (-l)"-iBo(n - l)!(5n2 - - 6^) /(m^ - (a^ + 6^))^ - 40^6^ 

m=l ^ 

20n + 5 35n^-35n3(a2-t-62)+4n(3a^ + 36^-4a262) 

^ 2(5n2 - 2a2 - 262) "+ 4(5n2 - 2a2 - 262) 

where is a solution of the second order difference equation 

4(4n + 3)(4n + 5) (5^2 - 2a^ - 2b^)Ln+i + 2{n + l)p{n)Ln 
- n{n + l)(5(n + 1)^ - 2a^ - 2b'^)q{n)Ln-i = 0, n = 1, 2, . . . 
with initial conditions Lq = Cq, 

= Q - ^(a^ + 6^)) Ao + Q(a^ + 6^) - l(a^ + 6^ - 4a^6^) - 1) Co, 

and 
(4) 

p(n) = 30n^ + 105n^ + n\U5 - 52(a^ + 6^)) + n\lOO - 130(0^ + 6^)) 
+ n^(35 - 124(a2 + 6^) + 56(a^ + 6^) - 2080^6^) + ^^(5 - 56(0^ + 6^) 
+ 84(a' + 6") - SUa^"") + n{80a%^{a'' + 6^) - 16(a^ + 6^) + 48(a^ + 6'* - 30^6^) 
- 14(a2 + 6^)) + (10(a2 - 6^)2 - 2(0^ + 6^) + 400^6^(02 + 6^) - 8(0^ + 6^)), 

q{n) = - 6^6(0' + 6=^) + n\9{a^ + b") + 30a%^) 

- n\28a^b\a^ + 6^) + 4(0^ + 6^)) + 16a%^{a^ - b^. 
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If in Theorem [1] we take Bq = 1, Aq = Cq = 0, then L„ = for all n > and putting 



,2 



X 



- 2 • - 2 

we get the identity ([3]). 

If Aq = 1, Bq = Cq = a = b = 0, we get the following series for ({4) mentioned by 
Markov in |10l p. 18]: 

n=l ^ 

where Lq = 0, Li = 1/3, and 

4(4n + 3)(4n + 5)L„+i + 2(n + 1)^(6^3 + + 5n + l)L„-n^(n + l)^L„_i = 0, n>l. 
Theorem 2. Let x,y be complex numbers such that |xp + |?/|^ < 1. Then 

k _ 1 ^ (-i)"-v(n) nr=\(K - a;')' + V) 



^A;4-x2P-^4 2^ n™= J^^ - ^^^^ - ' 

r (n) = 205n^ - IGOri^ + (32 - 62x^)n'^ + 40x^n^ + {x^ - 8x^ - 25y^)n^ + lOy^n + y^{x^ - 2) . 
Since 

fc=l ^ n=0 m=0 ^ ^ 

the formula ([7]) generates Apery-like series for all C{2n + 4m + 3), n, m > 0, convergent 
at the geometric rate with ratio 2^^°. So, for example, if x = ?/ = 0, we get Amdeberhan 
and Zeilberger's series [2] for ^(3) 

_ 1 >^ (-l)"-i(205n2 - 160n + 32) 

If y = 0, we recover Theorem 4 from [6]. If x = 0, we find, in particular, the following 
expression for C{7) : 

^n/'oc;^2 



..^,l^ i-m'^5n^-'^0n + 2) 

^ 2^ Ti9/'2n\5 



1 >^ (-l)"(205n^ - IGOn + 32) 1 ^\ 

n=l "' V n / \m=l m,=l 
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2. Proof of Theorem [H 
As usual, let (A)i/ be the Pochhammer symbol (or the shifted factorial) defined by 

^ r(A) |a(A + 1)...(A + z/-1), z/gN. 

Let a, b be complex numbers such that \a\ < 1, \b\ < 1. We start with the kernel 

(1 + a)fc(l - a)fc(l + - b)k 



(1 + a)„+/c+i(l — a)n+A;+l(l + — 

and define two functions 

F{n, k) = H{n, k){A{n) + B{n){k + 1) + C{n){k + 1)^), 
G{n, k) = H{n, k){D{n) + E{n)k + ir(n)A;2 + L{n)k^), 

with 7 unknown coefficients A{n), B{n), C{n), D{n), E{n), K{n), L{n) as functions of n. 
We require that 

(8) F{n + l,k)- F{n, k) = G{n, k + 1) ~ G{n, k). 

Since -F(n, k) and G{n, k) are not proper hypergeometric, the pair (F, G) is a Markov 
WZ-pair (see [TS^, p. 8], [T5] for definitions). 

Substituting F, G into ([H]) and cancelling common factors we get the following equa- 
tion of degree 6 in a variable k : 

{{n + k + 2f - a^){{n + k + 2f- b''){A{n) + B{n){k + 1) + G{n){k + 1)^) 

- A{n + 1) - B{n + l){k + 1) - C(r2 + l){k + 1)^ = ((n + A; + 2f - a^) 

X ((n + A; + 2f - b^){D{n) + E{n)k + K{n)e + L{n)k'') - {{k + if - a^) 

X ((A; + 1)2 - b^){D{n) + E{n){k + 1) + K{n){k + if + L{n){k + if). 

To satisfy condition ([H]), all the coefficients of the powers of {k + 1) in the equation 
(Q must vanish. This leads to a system of first order linear recurrence equations with 
polynomial coefficients for A{n), B{n), C{n), D{n), E{n), K{n), L{n) 

(10) G{n) = (4ni - 3)L(n), B{n) = (Aui - 2)K{n) - (lOnl - 3)L(n), ni = n + 1, 

(11) A{n) = (4^1 - l)E{n) - {lOnj - l)K{n) + (20n^ + 2ni{a^ + 6^) - l)L{n), 

(12) 4:D{n) = 10niE{n) - (20n? + 2a'^ + 2b^)K{n) + (35n? + llni{a^ + b^))L{n), 



2(4ni + l)L{n + 1) = 2ni(5n? - 20^ - 2b^)E{n) - 2n?(15n? - 6{a^ + b'^))K{n) 
^^^^ + ni(63nt - 17n?(a2 + 6^) _ 4(^4 + fe^^^^^^^^ 

2(4^1 + 2)K{n + 1) - 2(10n? + 20r2i + 7)L(n + 1) 
(14) = 2nl{hn\ - 2{a^ + b''))E{n) - 2{lQnl - Mlia" + 6^) + ^^(a^ - b'' f)K{n) 
+ (70^6 - 31n^(a2 + fe^) + ^2(3^4 ^ 3^4 _ 14^252))^^^)^ 
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4(20(^1 + 1)'^ + 2{ni + l)(a^ + b') - l)L{n + 1) - 4(10nf + 20ni + 9)K{n + 1) 
+ 4(4ni + 3)E{n + 1) = {6nl - 6nl{a^ + 6^) + 16a%^ni)E{n) 

^^^^ - {20nl - TlnXia" + \?) + 2n\{a^ + 6^ + 22a%^))K{n) 

+ {45nl - 48nl{a^ + 6^) + n?(3a^ + 36^ + 860^6^) + 8a%'^ni{a^ + b'^))L{n). 

Multiplying equation (fT3|) by ni and subtracting from we get 

2K{n + 1) - 7K + + 1) = - ^ 2(2n/+ l) " 

which yields 

/f(„) _ !„,i(„,, = (-l)"PA-(0)-7L(0))n! ^ /K - - - 4.^6^ 
^ ^ 2 ^ ^ 2"! J-l I 2m + 1 / 

m=l ^ ' 

From (Uni) it follows that 2fs:(0) = 5(0) + 7L(0) and therefore, 

(16) Kin) = -uMn) + ^ H ^^^^^^ J " 

m=l ^ ' 

Substituting (fT6ll into (fT3l) yields the formula for £'(n) 

4^1 + 1 42nf-25n2(a2 + 62)+4(a4 + 64) 
^(^) = — — 9 n T^Lin + 1) H ^ — -7^^ -Lin) 

^^^^ 35(0)(-l)"ni! A / (m2-a2- 62)2,4^2^2 

+ 11 ( 2m + 1 

m=l ^ 

Substituting f|T6l) and f|T71) into f|T2l) gives the formula for 

(18) 



_ (40r2i + 10)L(n + 1) + (35r2f - 35nf (a2 + 62) + 4ni(3a4 + 36'^ - 4a262))L(n) 
^""^ ~ 4(5^2 - 2a2 - 262) 

(-l)"5(0)n!(5n2 - a2 - 62) " /(m2 - a2 - 62)2 _ ^^2^2 
^ 2"i+i V 2m + 1 

m=l ^ 

Finally, substituting (fT6il . (fTTll into ( fTSl) gives the second-order difference equation for 
L{n) 

4(4n + 3)(4n + 5)(5n2 - 2a^ - 26^)L(ra + 1) + 2(n + \)'p{n)L{n) 
- n{n + l)(5(r2 + 1)^ - 20^ - 2\?)q{n)L{n - 1) = 0, 

with initial conditions i^(0) = C(0), 

L(l) = Q - A(a2 + ^2)^ ^(0) + + ^2) _ _|(^4 + ^4 _ _ 

derived from ( fTOl) . ( |TT1) . (|T3l) and polynomials p{n),q{n) defined in (jlj), ([5]). 
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If we put /(n) = L(n)/(n!)^, n = 0, 1, 2, . . . , then it is easily seen that the sequence 
l{n) satisfies the following recurrence equation: 

4(4n + 3)(4n + 5){5n^ - 2a^ - 2h^)n\n + lfl{n + 1) + 2r?p{n)l{n) 

- (5(n + 1)^ - 2a^ - 2l?)q{n)l{n - 1) = 0. 

Its characteristic polynomial 64A^ + 12A — 1 = has two different zeros Ai = —1/4, 
A2 = 1/16, and by Poincare's theorem, we get 

(18) ^^im)'<i 



n— »oo 



n\Y I ~ 4 



The limit inequality (]T9l) implies 



00 00 
lim V A;) = 0, limVG(n,fc) = 

n^oo ' ^ fc— >oo ' ^ 
k=0 n=0 

and therefore, we have (see |12J p. 9], [8, §2]) 

00 00 
Y,F{0,k) = Y,G{n,0), 

k=0 n=0 

yielding the theorem with dn = D{n — 1), L„ = L{n), Aq = A{0), Bq = B{0), Cq = 
C{0). □ 

3. Proof of Theorem [21 

To deduce (ED from take A{0) = C(0) = 0, B{0) = 1 and apply the following 
formula (see [HI Corollary 2]): 

00 00 

(20) Yl ^(0' ^) = ^) + ^('^' ^ + !))• 

fc=0 n=0 

Since in this case L{n) = for all n > 0, an easy computation of the right-hand side of 
([20]) by ([lOD, (HI]), ([T6])-([I8D and substitution lead to the desired conclusion. □ 
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